Screening and identification of potential MERS-CoV papain-like protease (PLpro) inhibitors; Steady-state kinetic and Molecular dynamic studies.

MERS-CoV belongs to the coronavirus group. Recent years have seen a rash of coronavirus epidemics. In June 2012, MERS-CoV was discovered in the Kingdom of Saudi Arabia, with 2,591 MERSA cases confirmed by lab tests by the end of August 2022 and 894 deaths at a case-fatality ratio (CFR) of 34.5% documented worldwide. Saudi Arabia reported the majority of these cases, with 2,184 cases and 813 deaths (CFR: 37.2%), necessitating a thorough understanding of the molecular machinery of MERS-CoV. To develop antiviral medicines, illustrative investigation of the protein in coronavirus subunits are required to increase our understanding of the subject. In this study, recombinant expression and purification of MERS-CoV (PLpro), a primary goal for the development of 22 new inhibitors, were completed using a high throughput screening methodology that employed fragment-based libraries in conjunction with structure-based virtual screening. Compounds 2, 7, and 20, showed significant biological activity. Moreover, a docking analysis revealed that the three compounds had favorable binding mood and binding free energy. Molecular dynamic simulation demonstrated the stability of compound 2 (2-((Benzimidazol-2-yl) thio)-1-arylethan-1-ones) the strongest inhibitory activity against the PLpro enzyme. In addition, disubstitutions at the meta and para locations are the only substitutions that may boost the inhibitory action against PLpro. Compound 2 was chosen as a MERS-CoV PLpro inhibitor after passing absorption, distribution, metabolism, and excretion studies; however, further investigations are required.

Multi-ligand molecular docking, simulation, free energy calculations and wavelet analysis of the synergistic effects between natural compounds baicalein and cubebin for the inhibition of the main protease of SARS-CoV-2.

Combination drugs have been used for several diseases for many years since they produce better therapeutic effects. However, it is still a challenge to discover candidates to form a combination drug. This study aimed to investigate whether using a comprehensive in silico approach to identify novel combination drugs from a Chinese herbal formula is an appropriate and creative strategy. We, therefore, used Toujie Quwen Granules for the main protease (Mpro) of SARS-CoV-2 as an example. We first used molecular docking to identify molecular components of the formula which may inhibit Mpro. Baicalein (HQA004) is the most favorable inhibitory ligand. We also identified a ligand from the other component, cubebin (CHA008), which may act to support the proposed HQA004 inhibitor. Molecular dynamics simulations were then performed to further elucidate the possible mechanism of inhibition by HQA004 and synergistic bioactivity conferred by CHA008. HQA004 bound strongly at the active site and that CHA008 enhanced the contacts between HQA004 and Mpro. However, CHA008 also dynamically interacted at multiple sites, and continued to enhance the stability of HQA004 despite diffusion to a distant site. We proposed that HQA004 acted as a possible inhibitor, and CHA008 served to enhance its effects via allosteric effects at two sites. Additionally, our novel wavelet analysis showed that as a result of CHA008 binding, the dynamics and structure of Mpro were observed to have more subtle changes, demonstrating that the inter-residue contacts within Mpro were disrupted by the synergistic ligand. This work highlighted the molecular mechanism of synergistic effects between different herbs as a result of allosteric crosstalk between two ligands at a protein target, as well as revealed that using the multi-ligand molecular docking, simulation, free energy calculations and wavelet analysis to discover novel combination drugs from a Chinese herbal remedy is an innovative pathway.

Glycyrrhizin as a promising kryptonite against SARS-CoV-2: Clinical, experimental, and theoretical evidences.

COVID-19 is the most devastating disease in recent times affecting most people globally. The higher rate of transmissibility and mutations of SARS-CoV-2 along with the lack of potential therapeutics has made it a global crisis. Potential molecules from natural sources could be a fruitful remedy to combat COVID-19. This systematic review highlights the detailed therapeutic implication of naturally occurring glycyrrhizin and its related derivatives against COVID-19. Glycyrrhizin has already been established for blocking different biomolecular targets related to the SARS-CoV-2 replication cycle. In this article, several experimental and theoretical evidences of glycyrrhizin and related derivatives have been discussed in detail to evaluate their potential as a promising therapeutic strategy against COVID-19. Moreover, the implication of glycyrrhizin in traditional Chinese medicines for alleviating the symptoms of COVID-19 has been reviewed. The potential role of glycyrrhizin and related compounds in affecting various stages of the SARS-CoV-2 life cycle has also been discussed in detail. Derivatization of glycyrrhizin for designing potential lead compounds along with combination therapy with other anti-SARS-CoV-2 agents followed by extensive evaluation may assist in the formulation of novel anti-coronaviral therapy for better treatment to combat COVID-19.

Interactive deciphering electron-shuttling characteristics of Coffea arabica leaves and potential bioenergy-steered anti-SARS-CoV-2 RdRp inhibitor via microbial fuel cells.

Due to the pandemics of COVID-19, herbal medicine has recently been explored for possible antiviral treatment and prevention via novel platform of microbial fuel cells. It was revealed that Coffea arabica leaves was very appropriate for anti-COVID-19 drug development. Antioxidant and anti-inflammatory tests exhibited the most promising activities for C. arabica ethanol extracts and drying approaches were implemented on the leaf samples prior to ethanol extraction. Ethanol extracts of C. arabica leaves were applied to bioenergy evaluation via DC-MFCs, clearly revealing that air-dried leaves (CA-A-EtOH) exhibited the highest bioenergy-stimulating capabilities (ca. 2.72 fold of power amplification to the blank). Furthermore, molecular docking analysis was implemented to decipher the potential of C. arabica leaves metabolites. Chlorogenic acid (-6.5 kcal/mol) owned the highest binding affinity with RdRp of SARS-CoV-2, showing a much lower average RMSF value than an apoprotein. This study suggested C. arabica leaves as an encouraging medicinal herb against SARS-CoV-2.

Identification of alkaloids from Justicia adhatoda as potent SARS CoV-2 main protease inhibitors: An in silico perspective.

The COVID-19 pandemic, caused by SARS CoV-2, is responsible for millions of death worldwide. No approved/proper therapeutics is currently available which can effectively combat this outbreak. Several attempts have been undertaken in the search of effective drugs to control the spread of SARS CoV-2 infection. The main protease (Mpro), key component for the cleavage of the viral polyprotein, is considered to be one of the important drug targets for treating COVID-19. Various phytochemicals, including polyphenols and alkaloids, have been proposed as potent inhibitors of Mpro. The alkaloids from leaf extracts of Justicia adhatoda have also been reported to possess anti-viral activity. But whether these alkaloids exhibit any inhibitory effect on SARS CoV-2 Mpro is far from clear. To explore this in detail, we have adopted computational approaches. Justicia adhatoda alkaloids possessing proper drug-likeness properties and two anti-HIV drugs (lopinavir and darunavir; having binding affinity -7.3 to -7.4 kcal/mol) were docked against SARS CoV-2 Mpro to study their binding properties. Only one alkaloid (anisotine) had interaction with both the catalytic residues (His41 and Cys145) of Mpro and exhibited good binding affinity (-7.9 kcal/mol). Molecular dynamic simulations (100 ns) revealed that Mpro-anisotine complex is more stable, conformationally less fluctuated; slightly less compact and marginally expanded than Mpro-darunavir/lopinavir complex. Even the number of intermolecular H-bonds and MM-GBSA analysis suggested that anisotine is a more potent Mpro inhibitor than the two previously recommended antiviral drugs (lopinavir and darunavir) and may evolve as a promising anti-COVID-19 drug if proven in animal experiments and on patients.

Benchmarked molecular docking integrated molecular dynamics stability analysis for prediction of SARS-CoV-2 papain-like protease inhibition by olive secoiridoids.

Objectives: We performed a virtual screening of olive secoiridoids of the OliveNetTM library to predict SARS-CoV-2 PLpro inhibition. Benchmarked molecular docking protocol that evaluated the performance of two docking programs was applied to execute virtual screening. Molecular dynamics stability analysis of the top-ranked olive secoiridoid docked to PLpro was also carried out. Methods: Benchmarking virtual screening used two freely available docking programs, AutoDock Vina 1.1.2. and AutoDock 4.2.1. for molecular docking of olive secoiridoids to a single PLpro structure. Screening also included benchmark structures of known active and decoy molecules from the DEKOIS 2.0 library. Based on the predicted binding energies, the docking programs ranked the screened molecules. We applied the usual performance evaluation metrices to evaluate the docking programs using the predicted ranks. Molecular dynamics of the top-ranked olive secoiridoid bound to PLpro and computation of MM-GBSA energy using three iterations during the last 50 ps of the analysis of the dynamics in Desmond supported the stability prediction. Results and discussions: Predictiveness curves suggested that AutoDock Vina has a better predictive ability than AutoDock, although there was a moderate correlation between the active molecules rankings (Kendall's correlation of rank (τ) = 0.581). Interestingly, two same molecules, Demethyloleuropein aglycone, and Oleuroside enriched the top 1 % ranked olive secoiridoids predicted by both programs. Demethyloleuropein aglycone bound to PLpro obtained by docking in AutoDock Vina when analyzed for stability by molecular dynamics simulation for 50 ns displayed an RMSD, RMSF<2 Å, and MM-GBSA energy of -94.54 ± 6.05 kcal/mol indicating good stability. Molecular dynamics also revealed the interactions of Demethyloleuropein aglycone with binding sites 2 and 3 of PLpro, suggesting a potent inhibition. In addition, for 98 % of the simulation time, two phenolic hydroxy groups of Demethyloleuropein aglycone maintained two hydrogen bonds with Asp302 of PLpro, specifying the significance of the groups in receptor binding. Conclusion: AutoDock Vina retrieved the active molecules accurately and predicted Demethyloleuropein aglycone as the best inhibitor of PLpro. The Arabian diet consisting of olive products rich in secoiridoids benefits from the PLpro inhibition property and reduces the risk of viral infection.

Elucidating the enhanced binding affinity of a double mutant SP-D with trimannose on the influenza A virus using molecular dynamics.

Surfactant protein D (SP-D) is an essential component of the human pulmonary surfactant system, which is crucial in the innate immune response against glycan-containing pathogens, including Influenza A viruses (IAV) and SARS-CoV-2. Previous studies have shown that wild-type (WT) SP-D can bind IAV but exhibits poor antiviral activities. However, a double mutant (DM) SP-D consisting of two point mutations (Asp325Ala and Arg343Val) inhibits IAV more potently. Presently, the structural mechanisms behind the point mutations' effects on SP-D's binding affinity with viral surface glycans are not fully understood. Here we use microsecond-scale, full-atomistic molecular dynamics (MD) simulations to understand the molecular mechanism of mutation-induced SP-D's higher antiviral activity. We find that the Asp325Ala mutation promotes a trimannose conformational change to a more stable state. Arg343Val increases the binding with trimannose by increasing the hydrogen bonding interaction with Glu333. Free energy perturbation (FEP) binding free energy calculations indicate that the Arg343Val mutation contributes more to the increase of SP-D's binding affinity with trimannose than Asp325Ala. This study provides a molecular-level exploration of how the two mutations increase SP-D binding affinity with trimannose, which is vital for further developing preventative strategies for related diseases.

Unsymmetrical aromatic disulfides as SARS-CoV-2 Mpro inhibitors: Molecular docking, molecular dynamics, and ADME scoring investigations.

COVID-19 pandemic caused by very severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) agent is an ongoing major global health concern. The disease has caused more than 452 million affected cases and more than 6 million death worldwide. Hence, there is an urgency to search for possible medications and drug treatments. There are no approved drugs available to treat COVID-19 yet, although several vaccine candidates are already available and some of them are listed for emergency use by the world health organization (WHO). Identifying a potential drug candidate may make a significant contribution to control the expansion of COVID-19. The in vitro biological activity of asymmetric disulfides against coronavirus through the inhibition of SARS-CoV-2 main protease (Mpro) protein was reported. Due to the lack of convincing evidence those asymmetric disulfides have favorable pharmacological properties for the clinical treatment of Coronavirus, in silico evaluation should be performed to assess the potential of these compounds to inhibit the SARS-CoV-2 Mpro. In this context, we report herein the molecular docking for a series of 40 unsymmetrical aromatic disulfides as SARS-CoV-2 Mpro inhibitor. The optimal binding features of disulfides within the binding pocket of SARS-CoV-2 endoribonuclease protein (Protein Data Bank [PDB]: 6LU7) was described. Studied compounds were ranked for potential effectiveness, and those have shown high molecular docking scores were proposed as novel drug candidates against SARS-CoV-2. Moreover, the outcomes of drug similarity and ADME (Absorption, Distribution, Metabolism, and Excretion) analyses have may have the effectiveness of acting as medicines, and would be of interest as promising starting point for designing compounds against SARS-CoV-2. Finally, the stability of these three compounds in the complex with Mpro was validated through molecular dynamics (MD) simulation, in which they displayed stable trajectory and molecular properties with a consistent interaction profile.

Montelukast is a dual-purpose inhibitor of SARS-CoV-2 infection and virus-induced IL-6 expression identified by structure-based drug repurposing.

Drug-repurposing has been instrumental to identify drugs preventing SARS-CoV-2 replication or attenuating the disease course of COVID-19. Here, we identify through structure-based drug-repurposing a dual-purpose inhibitor of SARS-CoV-2 infection and of IL-6 production by immune cells. We created a computational structure model of the receptor binding domain (RBD) of the SARS-CoV-2 spike 1 protein, and used this model for insilico screening against a library of 6171 molecularly defined binding-sites from drug molecules. Molecular dynamics simulation of candidate molecules with high RBD binding-scores in docking analysis predicted montelukast, an antagonist of the cysteinyl-leukotriene-receptor, to disturb the RBD structure, and infection experiments demonstrated inhibition of SARS-CoV-2 infection, although montelukast binding was outside the ACE2-binding site. Molecular dynamics simulation of SARS-CoV-2 variant RBDs correctly predicted interference of montelukast with infection by the beta but not the more infectious alpha variant. With distinct binding sites for RBD and the leukotriene receptor, montelukast also prevented SARS-CoV-2-induced IL-6 release from immune cells. The inhibition of SARS-CoV-2 infection through a molecule binding distal to the ACE-binding site of the RBD points towards an allosteric mechanism that is not conserved in the more infectious alpha and delta SARS-CoV-2 variants.

In search of RdRp and Mpro inhibitors against SARS CoV-2: Molecular docking, molecular dynamic simulations and ADMET analysis.

Corona Virus Disease 2019 (COVID-19) caused by Severe Acute Respiratory Syndrome coronavirus (SARS CoV-2) has been declared a worldwide pandemic by WHO recently. The complete understanding of the complex genomic structure of SARS CoV-2 has enabled the use of computational tools in search of SARS CoV-2 inhibitors against the multiple proteins responsible for its entry and multiplication in human cells. With this endeavor, 177 natural, anti-viral chemical entities and their derivatives, selected through the critical analysis of the literatures, were studied using pharmacophore screening followed by molecular docking against RNA dependent RNA polymerase and main protease. The identified hits have been subjected to molecular dynamic simulations to study the stability of ligand-protein complexes followed by ADMET analysis and Lipinski filters to confirm their drug likeliness. It has led to an important start point in the drug discovery and development of therapeutic agents against SARS CoV-2.

Screening of inhibitors against SARS-CoV-2 spike protein and their capability to block the viral entry mechanism: A viroinformatics study.

SARS-CoV-2, previously named 2019 novel coronavirus (2019-nCoV), has been associated with the global pandemic of acute respiratory distress syndrome. First reported in December 2019 in the Wuhan province of China, this new RNA virus has several folds higher transmission among humans than its other family member (SARS-CoV and MERS-CoV). The SARS-CoV-2 spike receptor-binding domain (RBD) is the region mediating the binding of the virus to host cells via Angiotensin-converting enzyme 2 (ACE2), a critical step of viral. Here in this study, we have utilized in silico approach for the virtual screening of antiviral library extracted from the Asinex database against the Receptor binding domain (RBD) of the S1 subunit of the SARS-CoV-2 spike glycoprotein. Further, the molecules were ranked based on their binding affinity against RBD, and the top 15 molecules were selected. The affinity of these selected molecules to interrupt the ACE2-Spike interaction was also studied. It was found that the chosen molecules were demonstrating excellent binding affinity against spike protein, and these molecules were also very effectively interrupting the ACE2-RBD interaction. Furthermore, molecular dynamics (MD) simulation studies were utilized to investigate the top 3 selected molecules' stability in the ACE2-RBD complexes. To the best of our knowledge, this is the first study where molecules' inhibitory potential against the Receptor binding domain (RBD) of the S1 subunit of the SARS-CoV-2 spike glycoprotein and their inhibitory potential against the ACE2-Spike has been studied. We believe that these compounds can be further tested as a potential therapeutic option against COVID-19.

In silico evaluation of flavonoids as effective antiviral agents on the spike glycoprotein of SARS-CoV-2.

The novel coronavirus pandemic has spread over in 213 countries as of July 2020. Approximately 12 million people have been infected so far according to the reports from World Health Organization (WHO). Preventive measures are being taken globally to avoid the rapid spread of virus. In the current study, an in silico approach is carried out as a means of inhibiting the spike protein of the novel coronavirus by flavonoids from natural sources that possess both antiviral and anti-inflammatory properties. The methodology is focused on molecular docking of 10 flavonoid compounds that are docked with the spike protein of SARS-CoV-2, to determine the highest binding affinity at the binding site. Molecular dynamics simulation was carried out with the flavonoid-protein complex showing the highest binding affinity and highest interactions. The flavonoid naringin showed the least binding energy of -9.8 Kcal/mol with the spike protein which was compared with the standard drug, dexamethasone which is being repurposed to treat critically ill patients. MD simulation was carried out on naringin-spike protein complex for their conformational stability in the active site of the novel coronavirus spike protein. The RMSD of the complex appeared to be more stable when compared to that of the protein from 0.2 nm to 0.4 nm. With the aid of this in silico approach further in vitro studies can be carried out on these flavonoids against the novel coronavirus as a means of viral protein inhibitors.